Semiconductor quantum dots in high magnetic fields: The composite-fermion view 
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We review and extend the composite fermion theory for semiconductor quantum dots in high 
magnetic fields. The mean-field model of composite fermions is unsatisfactory for the qualitative 
physics at high angular momenta. Extensive numerical calculations demonstrate that the micro- 
scopic CF theory, which incorporates interactions between composite fermions, provides an excellent 
qualitative and quantitative account of the quantum dot ground state down to the largest angular 
momenta studied, and allows systematic improvements by inclusion of mixing between composite 
fermion Landau levels (called A levels). 
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I. INTRODUCTION 

The system of interacting electrons confined 
to a two dimensional quantum dot and ex- 
posed to a strong magnetic field has been a 
subject of intense theoretical study for over two 

decades .l-2-3.4.5.6.7.8 .9.10.11.12.13.14.15.16.17.18.19.20.21.22.23.24.25 

Such quantum dots have been realized and studied in 
the laboratoryi2Si21iS^i2^ Exact diagonalization studies 
show that the ground states are strongly correlated, 
and the aim of theory is to achieve a satisfactory 
understanding of the correlations. It is also of interest to 
understand how this ties into our understanding of the 
FQHE,"^° obtained in the thermodynamic limit without 
confinement. 

The CF theory has been applied to parabolic quan- 
tum dots subjected to a strong magnetic field. The 
plot of ground state energy as a function of the angu- 
lar momentum (L) has a rich structure. In particular, 
downward cusps appear at certain values of L, which are 
consequently especially favorable. Early studies^'^ d^d^ 
demonstrated the CF theory to be promising. Specif- 
ically, a "mean-field model," in which the composite 
fermions are taken as noninteracting particles at an ef- 
fective angular momentum L* , with their mass or the cy- 
clotron energy treated as a phenomenological parameter, 
predicts cusps in the energy at certain magic angular mo- 
menta; these predictions are in agreement with the actual 
cusp positions in exact diagonalization studies at small 
angular momenta L, but discrepancies appear at large 
^ j^jis Puj-ther worls^ii^^ showed that these discrepancies 
are special to the mean-field model of the CF theory. A 
perfect agreement between the actual and the predicted 
cusp positions was obtained when the CF energies were 
calculated from microscopic wave functions. One of the 
surprising aspects was the success of the CF theory even 
at the largest angular momenta studied, which appears, 
at first, to be at odds with the classical crystal-like corre- 
lations found in exact diagonalization studies.— While 
both composite fermions and the crystal are generated 
by the repulsive interaction between electrons, the im- 



plicit assumption had been that one excluded the other. 
The work in Ref. ji^ showed that no logical inconsistency 
exists between the simultaneous formations of composite 
fermions and crystal-like structures at low fillings, and, 
furthermore, the formation of composite fermions itself 
induces crystal structure at low fillings. This crystal has 
been shown to be very well described as a crystal of com- 
posite fermions fS^i^ 

The aim of this paper is to review and extend the CF 
theory of quantum dots in high magnetic fields, and also 
provide many details left out in earlier papers. Section II 
briefly outlines the basics of the CF theory for quantum 
dot states, giving explicit wave functions for some simple 
cases. Section III describes the numerical methods (ex- 
act diagonalization, Lanczos, and CF diagonalization). 
The mean-field CF model is discussed in Sec. IV, and 
the "zeroth-order" CF diagonalization in Sec. V. Section 
VI illustrates how the results are improved by going to 
higher orders in the CF theory. The paper is concluded 
in Sec. VII. 



II. THE COMPOSITE-FERMION BASICS FOR 
QUANTUM DOTS 

Following the standard practice, we assume below 
parabolic confinement. This should be a good approx- 
imation for most quantum dots at low energies, and sim- 
plifies calculations because of the availability of exact so- 
lutions for single particle eigenstates. The Hamiltonian 
of interest is 
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where is the band mass of the electron, ujo is a mea- 
sure of the strength of the confinement, e is the dielectric 
constant of the host semiconductor, and rjk = \rj — r^j. 
We will specialize to the case of very large magnetic fields, 
when uJc = eB/mjjC S> ujq. Only the lowest Landau level 
(LL) is relevant in this limit. In that limit, at each an- 
gular momentum the eigenenergy neatly separates into a 
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trons in the lowest LL at angular momentum L is mapped 
into the noninteracting electron state at L* , where 



L* +pN{N - 1) , 



(3) 



N is the number of electrons, and p is an integer. The 
wave functions 



j<k 



(4) 



give ansatz wave functions for interacting electrons at L 
in terms of the known wave functions of noninteracting 
electrons at L*. Here are the wave functions for 
noninteracting electrons at L* (which in general occupy 
several Landau levels), a ~ 1,2,- •• , D* labels the dif- 
ferent states, Zj = Xj — iyj denotes the position of the 
jth electron, T^lll indicates projection into the lowest 
LL, are basis functions for interacting electrons at 
L, and D* is the dimension of the CF basis. We will 
restrict to states with the lowest kinetic energy at 
L*, and choose p so as to have the smallest dimension for 
the basis. The composite fermions carrying 2p vortices 
are labeled ^^CF's. At certain values of L, the above 
prescription produces only one state {D* = 1), which 
is the CF theory's answer for the ground state. In the 
notation of Ref. 10, this is a compact state, denoted by 
[iVo, A^i, • • • ], with Nj composite fermions compactly oc- 
cupying the innermost angular momentum orbitals of the 
jth CF level. At other values of L, when there are many 
CF basis states {D* > 1), we diagonalize the Coulomb 
Hamiltonian in the CF basis to obtain the ground state, 
using methods described earlierjiii^^ For any N, there 
are many values of L where the CF theory gives a unique 
answer, but in general, D* increases with TV. 



A. Examples of CF bases 



FIG. 1: Schematic depiction of (a) [N, 0], (b) [iV - 1, 1], and 
(c) [A^ — 2, 2] states in terms of composite fermions carrying 
two vortices (2p = 2). The composite fermions are shown 
as electrons with two arrows, where the arrows represent the 
vortices bound to composite fermions. The labels n and / 
denote the A level (CF Landau level) index and the angular 
momentum of the composite fermion, respectively. 



confinement part and an interaction part: 
E{L)^E,{L) + V{L), 



(2) 



where Ec{L) ~ {h/2)['[l — LUc]L, relative to the lowest LL, 
with = cj^ + 4:UJq, and V{L) is the interaction energy 
of electrons without confinement, but with the magnetic 
length replaced by an effective magnetic length given by 
£ = y^h/rribil. In the following, we will consider only V 
as a function of the angular momentum L; it must be 
remembered, however, that the confinement part must 
be added to determine the global ground state. 

In the CF theory^ i^°i^^ the interacting state of elec- 



Construction of the CF basis is in principle straight- 
forward. We consider some explicit examples. 

In many cases, the CF basis consists of a unique wave 
function, which simplifies the analysis tremendously. The 
simplest wave function is the one for one filled A level (CF 
Landau level), shown schematically in Fig. [TJa): 



A 



N 

n 

,i=l 

j<k 



2p 



Zk) 



2p+ip-2:"ikiiV4 



(5) 



where A denotes an antisymmetric Slater determinant, 
i.e. 
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It corresponds to L = (2p+l)N{N -l)/2. It is identical 
to the Laughlin wave function. 

The state [N - 1,1] at L = {2p+l)N{N - l)/2-N, 
shown in Fig. [lib) , represents the wave function 
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FIG. 2: Schematic depiction of basis states at (a) L — 33; (b) 
L = 34 for A'' = 6 particles. 



and 
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This state has been interpreted as a charged quasiparticle 
excitation of the = l/(2p+l) FQHE state,^ The [N - 
2,2] state (Fig. [He)) occurs at L = {2p+l)N{N-l)/2- 
2{N — 1), and has the wave function 
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For many angular momenta, the basis contains more 
than one state. For six composite fermions at L = 33 
we have two compact states [3,3] and [4,1,1], as depicted 
schematically in Fig. [^Ja). The basis functions are given 
by 



L=33 



,1-3 



i=3 
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These have the same CF kinetic energy, leading to D* = 
2 at L = 33. The five basis states at L = 34 are shown in 
Fig. [2Kb), derived from either [4,1,1] or [3,3] by increasing 
the angular momentum by one unit. The corresponding 
basis functions are written as 
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Here 



vUz) ^ z'Liizz*/2) 
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with LI^{x) being an associated Laguerre polynomial. 

The projection into the lowest LL is accomplished by 
the method outlined in the literature. To give a simple 



example, consider the state [N — 1,1] state in Eq. ((< 
We use the identity (apart from the constant factor), 
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TABLE I: Compaxison between the CF and the exact energies (Vcf and Vix) for iV = 4. The symbols in this and the following 
tables have the following meaning: L is the angular momentum; D is the dimension of the full lowest LL Fock space; and 
D* is the dimension of the CF basis. The energies are in units of e^/e£, where e is the dielectric constant of the background 
semiconductor and £ = -^/hc/eB is the magnetic length. 
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with 



Ji = [| {Zi - Zj] 



and project each element to the lowest LL, resulting in 
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(14) (15) 
The final step is achieved by placing all of the z* 's to the 
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TABLE III: Comparison between the CP and the exact energies (Vcf and Vex) for TV = 6. 



left, and substituting 2d/dzj into z* with the assump- 
tion that the derivatives do not act on the Gaussian fac- 
^Qj, ^35,38 rpj^g resulting wave function is then given by 



^ = A 



1 



N 

j^l 3 1=2 



,1-2 



(16) 



Lowest LL projection of other can be accomplished sim- 
ilarly, although the details are more complicated. 



routines for small L, or a modified Lanczos algorithm 
for larger L. In either cases, it is essential to know the 
Coulomb matrix elements. In the second quantization 
language, the Coulomb Hamiltonian is written as 



2 ^ 



(r, s\V\t, u) a], al a„ a*, 



(17) 



III. NUMERICAL METHODS 

A. Exact diagonalization 

We calculate the exact interaction energy for small sys- 
tems by either numerical diagonalization using standard 



where the operator al. {ar) creates (annihilates) an par- 
ticle at state \r) with angular momentum r. V is the 
Coulomb interaction in real space and {r,s\V\t,u) is 
its restriction to the lowest Landau level Hilbert space. 
Several expressions for {r, s\V\t,u) exist in the litera- 
turei^i^'^ In this work we use the formula derived by 
TsiperiS: 



where 



{s + r, t\V\s., t + r)^ ——^-^^[A:, + B^, A^ (18) 

i^vO (r + .)!r(i+r + t + z)' ir + ^)W{|+r + t + ^)^2+' + ^'^■ ^^'^ 

I 



Each term in the sum on the right hand side of Eq. pop is 
positive definite, which makes this expression more stable 



in numerical calculations, which is especially important 
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TABLE IV: Comparison between the CF and the exact energies (Vcf and Vex) for N = 7. 
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TABLE V; Comparison between the CF and the exact energies (Vcf and Vox) for N — 8. 



for large L. 

Exact numerical diagonalization is limited to systems 
with small numbers of particles at small angular momen- 
tum L, because the size of the Hilbert space grows ex- 
ponentially fast with TV and L. For example, for N = 6, 
the Fock space dimension grows from 21 to 701661 as L 
is increased from 21 to 148. When D becomes large, we 
have obtained the ground state energy using the modi- 
fied Lanczos algorithm of Gagliano et al.^ Briefiy, the 
algorithm begins with an initial guess of the ground state 
I'f/'o)- By applying the Hamiltonian TL onto a new 

state I -01) is defined as the following 
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where the notation (7i") represents the expectation value 
('0o|'^"|^o)- It is straightforward to verify that 1-01) is 
normalized and orthogonal to |0o) by this construction. 
The Hamiltonian TL is now diagonalized in the subspace 
spanned by {|0o), IV'i)}; ^ind the lowest eigenvalue e and 
its corresponding eigenvector |0;i) are chosen as a better 
approximation the to true ground state and its energy 
than {H) and |0o)- In terms of (H") and li/'o); the "bet- 
ter" ground state energy and corresponding state can be 
written as 
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TABLE VI: Comparison between the CF and the exact energies (Vcf and Vex) for = 9. 
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TABLE VIL Comparison between the CF and the exact energies (Vcf and 14x) for N — 10. 



in which 

a - / - ^W+l, (23) 

^ 2{{n^) - (H)2)3/2 ' ^ > 

b = V(^^> - in)'- (25) 

The state jV'o) is then used as a new initial guess and 
the entire procedure is iterated until the relative energy 
difference |e— {H) \ is smaller than a predefined tolerance. 

We place the program on a node of a Beowulf type PC 
cluster, each node consists a dual Pentiumlll IGHz pro- 
cessor. The modified Lanczos algorithm converges rel- 
atively fast to the true ground state. The most time 
consuming part of the calculation is the Coulomb matrix 
element, especially when L is very large. For example, 
when iV = 6 it takes approximately 24 CPU hours to 
obtain the ground state energy for L — 135. The largest 
system we have studied by the Lanczos method has a 
Fock space dimension of 817789. We note that for large 
L, our energies are slightly lower than those in Ref. llSl . 
presumably because they work with a truncated basis. 



B. CF diagonalization 

CF diagonalization refers to a diagonalization of the 
full Coulomb Hamiltonian in a correlated CF basis. 
The basis wave functions correspond to states with low 
"CF kinetic energies" (the kinetic energy is given by 
^ '^jLi nj) with the restriction that the total angular 
momentum is L, i.e., 

N 

Y,lj+PN{N-1)=L, (26) 
i=i 

where Uj and Ij are the A-level and the angular mo- 
mentum indices, respectively, for composite fermions. In 
most of this paper, we choose states that have the lowest 
CF kinetic energy (zeroth order CF diagonalization). By 
allowing hybridization with CF states with higher kinetic 
energies (higher order CF diagonalization), more accu- 
rate approximations can be obtained, as seen explicitly 
below. As illustrated in the previous section, the wave 
function corresponding to the ath basis at L is then 



8 




FIG. 3: The exact interaction energy (circle) vs. the mean- 
field CF ground-state energy (square) as a function of the 

angular momentum L for N = Q. The exact energy is plotted 
in units of e^/e£ while the CF energy is in arbitrary units. 



given by 



V'„w4«)(^i)-^„(=),4=)(22) 

(a = l,2,...,Z)*) 



"JV ''JV 



where 



(27) 



(28) 



The CF basis states constructed in this manner are 

not orthogonal to one another at given i, and sometimes 
not even hnearly independent. We use Gram-Schmidt 




20 40 60 80 100 120 140 



(a) 




FIG. 4: The relative errors of CF ground-state energy com- 
pared with the exact energies for (a) N = 4; (b) iV = 5; (c) 
N = 6; (d) N = 7. 



procedure to orthogonahze the states. The orthogonal 
basis states |^) are expressed as 



l^a) = \Va) 



(29) 
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FIG. 5: The exact interaction energies V (dashes) are given as a function of the angular momentum L for (a) = 4 (b) 
= 6 particles. The dots are the predictions of the zeroth-order CF theory. Different panels correspond to L regions where 
composite fermions of different flavors are relevant. The energies are quoted in units of e^/e£. 



where the normalized state \ria) is defined by 



(30) 



From the relation in Eq. we can find the recursion 
relation for the transformation matrix Uaf3, defined by 



Ual3 



where O 



Q/3 



E 



1 





■Ujf3 for P < a, 

for P — a, 
for l3 > a, 



(31) 

iValvp)- The computation of Ua/3 en- 
ables us to calculate the Coulomb Hamiltonian matrix 
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elements Va0 in orthonormal basis sets, 



where 

(33) 
(34) 



7,(5 



The Coulomb Hamiltonian matrix in Eq. (j32p is diago- 
nalized to obtain the energy and the wave function of the 
ground state in the CF theory. 

The matrix elements Oa/s and {rjalVlrjj^) have been 
evaluated by the Metropolis Monte Carlo method. We 
have performed in access of 4 x 10^ Monte Carlo steps for 
each energy, which correctly gives five significant digits 
for < 7 and four significant digits for > 8. The 
error bar denotes the standard deviation obtained from 
four independent runs. 



IV. MEAN-FIELD MODEL 

It is customary to first consider the mean-field version 
of CF model, in which the interaction between composite 
fermions is assumed to vanish. Then the interaction en- 
ergy of electron systems is completely transformed into 
the "kinetic energy" of composite fermions. and can be 
evaluated in units of CF cyclotron energy huj* by sum- 
ming the A-level indices of all occupied CF states 



VcF- 



MF 



N 



huj* 



(35) 



Figure [3] compares the interaction energies predicted 
by the mean-field CF theory with the exact ground-state 
energies as a function of total angular momentum L for 
N — 6. For small angular momenta the mean-field CF 
model predicts correctly the positions of cusps on ground- 
state energy {e.g. see Fig. [3l^a)). However, with increas- 
ing L, it gradually fails to capture the correct positions of 
cusp states. Such discrepancies between the exact results 
and the mean-field predictions motivate the need for an 
investigation of the residual interactions between com- 
posite fermions, considered in the subsequent sections. 



V. ZEROTH-ORDER COMPOSITE FERMION 
THEORY 

Based on the CF diagonalization procedure described 
in Sec, nil Bl we have computed the CF ground state ener- 
gies at the zeroth order, i.e. by including only the lowest 
kinetic energy CF basis states. We have also carried out 
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FIG. 6: Exact energy spectrum for = 6. 



exact diagonalization for N < 10. It is possible to go to 
arbitrarily large values of L within the CF theory, but 
available computer memory and execution time restricts 
our exact diagonalization study to L < 148, 127, 118, 116, 
and 113 for TV = 6, 7, 8, 9, and 10, respectively. (To our 
knowledge ^ 34,42 ,43 ^j^g largest systems for which exact re- 
sults have been computed are N = 8, L = 140 with Z3ox = 
3023010 and TV = 10, L = 135 with D^^ = 2977866.) The 
results are listed in Tables Ul I VIII which show that the CF 
energies compare well to the exact ones. (Many of the 
exact energies have been reported previously in the liter- 
aturefiiii^i^^i^ but are included here for completeness.) 
The deviation between the CF and the exact energies 
(VcF and Vex) is L dependent but small in the entire 
range of L studied in this work. The largest deviations 
are 1.1% (L = 36,TV = 4), 0.6% {L = 102,7V = 5), 
0.5% (L = 135,7V = 6), 0.7% (L ^ 112, iV = 7), 0.9% 
{L = 89,N ^ 8), 1.3% {L = 114, iV = 9), and 0.2% 
{L ~ 80,N — 10). For = 4, the maximum deviation 
occurs for composite fermions carrying 6-8 vortices; for 
N > 5, the deviation grows generally with L in the range 
considered in this work (see Fig. 2]). 

The plot of interaction energy as a function of L in 
Fig. [5] gives a demonstration of the accuracy of the CF 
predictions. In addition to the quantitative accuracy, the 
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FIG. 7: Exact energies (dashes) and the CF energies (dots) for six particles as a function of the angular momentum L. The 
left, center, and right panels correspond to ^CFs, *CFs, and ®CFs, respectively. The CF energies are obtained in the zeroth 
order calculation (explained in text). For L > 62, the exact spectrum is truncated because the energies are obtained by the 
Lanczos method. The energies are multiplied by \/T, for convenience of illustration. 




20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 105 

L L L 

FIG. 8: Same as in Fig. [T] but with the CF energies obtained in the first order calculation. 



qualitative features of the energy versus L plot are repro- 
duced faithfully by the CF theory. For TV = 4, the major 
cusps occurs at L = 4n-|-2. The system with six particles 
exhibits more complicated features. For system of ^CFs 
(top panel) shows cusps at L = 21, 25, 27, 30, 33, 35, 39, 40 
and 45. As the flavor increases, the cusps at angular mo- 
menta other than L = bn become less prominent and 
eventually disappear. Such periodic behavior is consis- 
tent with the geometric interpretation] -'^^i^^ 



VI. NEXT ORDER CF THEORY 

The CF theory allows a systematic perturbative way 
of improving the results. Above we considered only the 
CF states with the lowest CF kinetic energy, called the 
zeroth-order CF theory. The next step is to include states 
with one more unit of the kinetic energy in the CF ba- 
sis^i The degree of improvement can be seen in Figs. [H 
[71 and [51 Figure [HI shows the full spectrum for six elec- 
trons in a range of angular momentum; Figure [3 shows 
the spectrum from the zeroth order CF diagonalization. 
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and Fig. [8] from the first order. Both the zeroth and 
first order spectra capture the quahtative behavior and 
the positions of the cusps, but the first order theory is 
quantitatively much more accurate. In both cases, the 
discrepancy between the CF and the exact energies (Vcf 
and Vex) grows with L, because while the dimension of 
the CF basis remains the same when L is changed by 
N{N — 1), the dimension of the full lowest LL Fock space 
increases rapidly. In spite of the small basis (which some- 
times contains only one state at the zeroth order), the CF 
theory is quantitatively satisfactory. 

VII. CONCLUDING REMARKS 

Several author a^^'^^ have noted that the CF theory 
fails to produce the cusp positions for large L. These 
comparisons, however, refer only to the mean-field model 
of composite fermions, in which the interaction energy 
of electrons at L is viewed as the kinetic energy of free 
fermions at an effective angular momentum L*, with the 
cyclotron energy treated as a parameter. The validity of 
the CF theory should not be confused with the validity of 
the mean-field picture, which serves, at best, as a useful 
guide; given its crudeness, it is in fact surprising it works 
as well as it does. A more substantive, microscopic test 
of the CF theory requires working with the correlated 
wave functions produced by the CF theory. 

Our extensive study of quantum dot states shows 
that the microscopic composite fermion theory, defined 
through wave functions, gives an excellent description in 
regions including both liquid-like and crystal-like ground 
states, and continues to be satisfactory from very low an- 
gular momenta to the largest angular momenta studied 
to date. It provides an accurate approximation for the 
ground state wave function and the ground state energy 
at every single L in the wide range studied, and correctly 
reproduces all cusps in plot of the ground state energy 
vs. L. Taken together, these results constitute a detailed 



verification for the validity of the composite fermion the- 
ory for quantum dots. 

It is expected, from general considerations, that the 
ground state at large L will resemble a classical crys- 
tal, because large L implies small density (or small filling 
factor) with particles far from one another, as a result of 
which the system behaves more or less classically. Ref- 
erence |1^ has studied a Hartree-Fock crystal trial wave 
function based on an analogy to the classical crystal 
ground state in a quantum dot. No crystalline corre- 
lations are put in by hand in our calculations described 
above, however. As implied by the successful compar- 
isons with the exact results, and also confirmed by an 
explicit calculation of the pair correlation functio n^^'^^ , 
the CF theory automatically generates a crystal of cor- 
rect symmetry. The CF approach offers many other ad- 
vantages over the Hartree-Fock electron crystal descrip- 
tion. The latter obtains wave functions and energies only 
for certain special values of L, and even then only for the 
ground state. The CF theory, on the other hand, pro- 
vides a quantitative understanding of states at all L. For 

= 6, Reference [l^ explicitly quotes energies from their 
approach for seven values of L in the range 75 < L < 135. 
For these angular momenta, the zeroth order CF theory 
gives lower energy in every case except at L = 135. (The 
state ai L = 135 has one filled A level, described by 
the Laughlin wave function.) At the first order, the CF 
results improve substantially for large L. As shown in 
Ref. [3^ an almost exact description of the crystallite at 
large L is obtained in terms of a crystal of composite 
fermions, wherein a combination of the crystal and CF 
physics are introduced right at the outset. It is gratifying 
that the principle that applies to the fractional quantum 
Hall effect in a bulk two-dimensional electron system also 
produces an understanding of the quantum dot physics 
at high magnetic fields. 
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